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We report on a study of the spin-1 ferromagnetic Bose-Einstein condensate with magnetic dipole- 
dipole interactions. By solving the non-local Gross-Pitaevskii equations for this system, we find three 
ground-state phases. Moreover, we show that a substantial orbital angular momentum accompanied 
by chiral symmetry breaking emerges spontaneously in a certain parameter regime. We predict that 
all these phases can be observed in the spin-1 *^Rb condensate by changing the number of atoms 
or the trap frequency. 
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The magnetic dipole-dipole interaction in ferromagnets 
is responsible for a rich variety of spin structures [l[ . Sim- 
ilar spin textures can also be expected to occur, due to 
the dipolar interaction, in ferromagnetic Bose-Einstein 
condensates (BECs), such as the spin-1 ^^Rb BEC. How- 
ever, spinor BECs differ from ferromagnets in that they 
exhibit spin-gauge symmetry that can generate mass flow 
by developing spin textures '2| . One may therefore won- 
der whether the dipole-induced spin texture can yield 
spontaneous mass current in the ground state. In this 
Letter we show that this is the case. 

An additional motivation for our work is the recent ob- 
servation of a dipolar BEC in a system of spin-polarized 
^^Cr atoms 3]. The ground state of spinor dipolar BECs 
has been studied by several researchers using a single- 
mode approximation 0, [^, 0]. Magnetism of dipolar 
BECs in one and two dimensional optical lattices has 
also been investigated 0, S]- In this Letter, we study the 
ground-state spin textures of a spinor dipolar BEC with- 
out invoking the single-mode approximation. We show 
the existence of three ground-state phases at zero mag- 
netic field. In particular, we identify chiral spin-vortex 
phase and show that spontaneous circulation with broken 
chiral symmetry emerges in this phase. 

We consider a system of N spin-1 atoms with mass 
M confined in a spin- independent potential Utra.p{r) = 
Muj'^{x^ + y-^ -\- z^)/2. The Hamiltonian of the system is 
given by [9| 

-fftot ^ J [■0j;,^(r)i?o^Ami(T-) 

+ Hdd, (1) 

where iprnif) is the annihilation operator of an atom 
in the magnetic sublevel m = 0, ±1 at point r, Hq = 



— ^^V^/(2M) + t/tiap(^), Fmm' represents the spin-1 
matrix, and ifdd is the contribution of the magnetic 
dipole-dipole interaction. The spin-independent and 
spin-dependent interactions are characterized by go — 
ATTh^{ao -f 2a2)/(3M) and gi = 4Trh'^{a2 - ao)/(3M), re- 
spectively, with fls (s = 0, 2) being the 5-wave scatter- 
ing length for the scattering channel with total spin s. 
For spin-1 *^Rb atoms we have < — gi ^ go, and the 
ground state is ferromagnetic. 

The magnetic dipole moment of an atom is given by 
p, = —gslJ^BS + giiJLT<fI, where /ib is the Bohr magne- 
ton, /.tN is the nuclear magneton, S and / are the elec- 
tronic and nuclear spin angular momenta, respectively, 
and gs and gi are the Lande g factors. The matrix ele- 
ment of the dipole moment between magnetic sublevels in 
the spin-1 hyperfine manifold is shown to be {m\jj\'m') — 
gplJ-BFmm' ■ For the case of ^^Rb with 5=1/2 and 
/ — 3/2, we have gp = {gs + ^gjfi^ / f^B)/'i — 1/2. Then 
the contribution of the magnetic dipole-dipole interaction 
to the Hamiltonian is given by 

^ i 7n ' ^ ) 

\r — r'\-^ 

where Cdd = Mo5|-Mb/(4''') with fio being the magnetic 
permeability of the vacuum, and e = (r — r')/\r — r'\. 
We ignore the coupling of the F = 1 manifold with the 
F = 2 manifold due to the dipolar interaction because 
the hyperfine splitting ^ 100 niK is much larger than the 
dipolar energy ~ 1 nK. 

Our system has three characteristic length scales: 
dipole healing length = h/^/2Mcddno, spin heal- 
ing length S^sp = fi'/ \/^^M\gi\no, and Thomas-Fermi 
(TF) radius Rtf — '2aho[go'no/{fiijj)]^^^ , where uq = 
huj/go[5ia2 -\- 2ao)7V/(8aho)]^/^ is the TF peak density 
and Oho = [h/{2Mu>)]^^'^. As shown below, the phase di- 
agram is characterized by two dimensionless parameters 
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Rrp/Ud and Rtf/£,sp- 

We investigate the mean-field ground state of the total 
Hamiltonian TJtot ■ The general form of the ferromagnetic 
spinor order parameter is given by 



e-*"cos2(/3/2) \ 
''^ y2sin(/3/2)cos(/3/2) , (3) 
e*"sin2(/3/2) / 



where n is the number density, 7 is the gauge, and 
s = (sin /3 cos a, sin /3 sin a, cos /3) is the unit vector par- 
allel to the magnetization. The dipolar interaction in- 
duces a spin texture characterized by a spatial depen- 
dence of Euler angles a, /?, 7. This dependence in- 
duces mass and spin velocity fields which can be cal- 
culated from continuity equations dn/dt + V • {nv) = 
and d{nSfj)/dt + V • [nv^^ + ns^^v) — 0, a.s v — 
~h/M{V-/ + cosPVa), = -h/{2M)s\v? l3Va, and 
<p + ivlp = -h/{2M)e'" [- sin/3cos/3Va -I- iV/3]. Here 
nv and nv^^ (/i = x, y, z) describe the current density 
of the mass and that of the /Lt-component of the spin 
vector s, respectively. In addition to the kinetic en- 
ergy Kq = J dr?i^/(2M)|V^/np arising from the spa- 
tial dependence of the total density, the system has 
kinetic energies i^mass — J rfrinM|Dp and i^spin = 
JdrnMj:^ Kp|2 = /drftV(4M)E^ jVs^P caused by 
the mass and spin currents, respectively. In the case of 
ferromagnets, the domain structure is mainly determined 
by the interplay between the dipolar interaction and the 
spin stiffness that has the form of i^spin- In the case of 
dipolar BECs, the additional term -fCmass is significant. 

The ground-state phases can be classified by the sym- 
metry of the order parameter. The total Hamilto- 
nian -fftot is invariant under space inversion Pi/'m(r) = 
4'm{—''') and time reversal T'0„j(r) — (— l)"'7/'l„(r). It 
also possesses rotational [i.e., S0(3)] symmetry in com- 
bined spin and coordinate space around an arbitrary 
axis. However, as we show below, the S0(3) symme- 
try is reduced to uniaxial symmetry and the obtained 
phases are eigenstates of the projected total angular mo- 
mentum Jz — —i{d/dip) + Fz in cylindrical coordinates 
r = (r, z), characterized by 



a{r) = ip + a{r, z), 
/3(r)=/?(r,z)=/3(r,-z), 
7(r) = -Ja{r) -|- 7(r, z). 



(4a) 
(4b) 
(4c) 



with J being the eigenvalue of Jz- We thus classify the 
obtained phases by the value of J, space-inversion sym- 
metry, and time-reversal symmetry. 

The ground-state phases are obtained by solving the 



non-local Gross-Pitaevskii (GP) equations, 
{Ho + gon)ipo + (51/+ - '^dd^+)-^ 



+ {gif- ~ "^dd^-)^^ = M^o, (5a) 



(iJo + 5o?^)^±l + (5i/t - Cddfo^) 



^0 
^/2 



±(91/2 - CA&hz)i;±\ = /iV'iij (5b) 

where [i is the chemical potential, n — IV'mP is the 
number density, and / = {fxify,fz) is the spin density 
defined by fz = IV'iP - IV'-iP and + ify = /+ = 
/* = \/2{iplipo + ipQip-i). When the system is purely 
ferromagnetic, the spin density satisfies / = ns. We 
define the scaled dipole field b — {bx,by,bz) as 6p(r) = 
-J2,.=x,y,zl dr'iSf,^ - 3e^,e^)U{r')/\r - r'\^ and b± = 
bx ± iby. The effective magnetic field produced by the 
surrounding magnetic dipoles is given by Cddb/(3_FA'B)- 

We numerically solve Eq. ^ in three dimensions by 
the imaginary-time propagation method. We prepare the 
initial state by first solving Eq. ([5]) with Cdd — and then 
randomizing the local spin directions while keeping the 
number density unchanged. Comparing the total energies 
of the obtained textures, we obtain the phase diagram in 
the parameter space of (i?TF/'Csp)^TF/Cdd) as shown in 
Fig. 1(a). For the spin-1 *^Rb BEG, we have 6p/Cdd = 
\/cdd/\gi\ = 0.30, and the system follows the dotted fine 
in Fig. 1(a). A typical spin configuration in each phase is 
shown in Figs. l(b)-l(d). The symmetry properties are 
summarized in Fig. 2. We now describe the properties of 
each phase. 

Polar-core vortex (PCV) phase. — This phase has to- 
tal angular momentum J = and has a polar-core, i.e., 
(V'i,'0o,V'-i) = 7^(0,1,0) at r = 0. The length of 
the spin vector \f\/n is less than unity for r < ^sp, and 
the spin configuration far from the core is described by 
a = ±7r/2, P = 7r/2, and 7 = const. Since the m = 1 
component has the phase winding number —1 and the 
TO = — 1 component has the phase winding number -1-1, 
the net mass current vanishes while the net spin current 
remains, giving a nonzero ifspin- Time-reversal symme- 
try and space-inversion symmetry are both broken, but 
the order parameter remains invariant under the com- 
bined PT operation. 

Flower (FL) phase 1T|. — This phase has J = 1 and can 
continuously develop from a spin-polarized state. The 
spin configuration is given by a{z > 0) = 0, a{z < 0) = 
TT, < /3 <C 7r/2, and 7 = const. Since /3 ~ 0, most of the 
total angular momentum Nh resides in the spin and the 
remaining small angular momentum is carried by mass 
current. The spin current is also small. The system has 
space-inversion symmetry while time-reversal symmetry 
is broken. 

Chiral spin-vortex (CSV) phase. — This phase has J — 
1 and the spin structure is described by < a < tt. 
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FIG. 1: (Color) (a) Phase diagram of a ferromagnetic dipolar 
BEC. The solid curve shows the first-order phase boundary 
between the J = and J = 1 phases and the broken line 
represents the second-order one and divides the phases with 
and without spin chirality, where PCV, CSV, and FL stand 
for polar-core vortex, chiral spin-vortex, and fiower phases, 
respectively. The schematic shown in each region represents 
the spin configuration (black arrows) and mass fiow (red ar- 
rows). The spin-1 ^'^Rb BEC (Csp/^dd = 0.30) traces the 
dotted line, which is referred to as the Rb line, as a function 
of iiTp/Csp- (b)-(d) Typical spin configuration in each phase. 
The top and bottom panels show the unit vector f{j')/\f{r)\ 
projected onto the x-z and x-y planes, respectively. The 
color of the arrows represents the magnitude of the spin den- 
sity normalized by the number density \f{r)\/n(r) (see the 
bottom scale). 
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FIG. 2: Broken (x) and unbroken (Q) symmetries for each 
phase in Fig. 1. The schematic in each cell shows the spin 
configuration after the transformation by P, T, and PT. 
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FIG. 3: Orbital angular momentum in each state as a function 
of i?TF/Cdd. The inset shows the ratio of the dipolar energy 
to the total energy E^d/Etot = {Hdd} / (Htot) for each state at 
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a{r, z) + 5(r, — z) ~ tt, and < P < tt. The gauge 7 also 
varies in space. As shown in Fig. 1(d), this phase has 
a coreless vortex whose size is of the order of ^dd- The 
system has substantial spin and mass currents around 
the rotational axis and has nonzero /Cmass and iCspin- 
Both time-reversal symmetry and space-inversion sym- 
metry are broken. 

The dependence of the dipolar energy is shown in the 
inset of Fig.O The dipolar energy is minimized when the 



spin texture has a flux-closure structure I4I, i.e., when 
the divergence of the magnetization is everywhere zero. 
The PCV state, which satisfies V • / = in all space, has 
the lowest dipolar energy of the three states. 

The CSV state is different from the FL state in that 
the chiral symmetry of the spin vortex is broken to form 
a flux-closure structure. When the CSV state has a 
nonzero total angular momentum in the +z direction, 
there are two degenerate whirling patterns of spins in 
the x-y plane: clockwise (5 < 0) and anticlockwise 
(a > 0). These patterns have the spin configurations 
of right-handed and left-handed screws, which are trans- 
formed into each other under space inversion P. The 
system selects one of the degenerate states by break- 
ing the chiral symmetry. Accompanied by the chiral 
symmetry breaking, the net orbital angular momentum 
Lz — J drJ2m^mi~^^9/^'fi)'^rn incrcascs drastically as 
shown in Fig. [3l where in each state is plotted as a 
function of i?TF/'^dd- The FL and CSV states both have 
a finite orbital angular momentum, but the dependence 
on i?TF/Cdd is quite different. Since the order parameters 
in the FL and CSV states are purely ferromagnetic, the 
spin texture depends only on i?TF /Cdd , and therefore Lz 
is independent of ^sp- 

Finally, we discuss possible experimental situations. 
To form a spin texture in the ground state, the Zeeman 
energy must be smaller than the dipolar interaction en- 



4 



0.5 p 

0.4 

0.3 

0.2 

0.1 







FL 




PCV 



1000 

( (0 [rad/s] N' 



2000 



from the Ministry of Education, Culture, Sports, Science 
and Technology of Japan. YK acknowledges support by 
the Japan Society for Promotion of Science (Project No. 
185451). MU acknowledges support by a CREST pro- 
gram of the JST. 



Note added. — The ground-state phase diagram of 
a spinor dipolar system has recently been discussed in 
Ref. [11. 



FIG. 4: Orbital angular momentum as a function of (luN^)^^^ 
for *^Rb along the Rb line in Fig. 1. 

ergy, otherwise the ground state is spin-polarized. This 
implies that B < 10'^ G for the spin-1 ^^Rb BEC with 
n = 10^^ cm"'^. The critical field increases in propor- 
tion to the atomic density. By changing the number of 
atoms or the trap frequency, all three phases can be ex- 
perimentally realized. Figure 2] shows the orbital angular 
momentum as a function of {cjN^y/^ along the Rb line 
in Fig. 1(a). In the case of a BEC with iV = 10*^, for ex- 
ample, the CSV phase exists in the region of 2tt x 70 Hz 
< a; < 27r X 630 Hz, where increases up to OANh. It is 
interesting to note that we can induce a phase transition 
by controlling the trap frequency. 

In summary, we have investigated the phase dia- 
gram of a spin-1 ferromagnetic Bose- Einstein condensate 
with magnetic dipole-dipole interactions and found three 
phases, namely, the polar-core vortex phase, the flower 
phase, and the chiral spin- vortex phases. The chiral spin- 
vortex phase has chirality in the formation of the spin 
vortex, and the topological spin structure spontaneously 
yields a substantial net orbital angular momentum. 
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